Coulomb Drag and High Resistivity Behavior in Double Layer 
Graphene 
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PACS 81 . 05 . lie - Graphene 

PACS 72 . 80 . Vp - Electronic transport in graphene 

Abstract. - We show that Coulomb drag in ultra-clean graphene double layers can be used for 
controlling the on/off ratio for current flow by tunning the external gate voltage. Hence, although 
graphene remains semi-metallic, the double layer graphene system can be tuned from conductive 
to a highly resistive state. We show that our results explain previous data of Coulomb drag in 
double layer graphene samples in disordered Si02 substrates. 



Introduction. — Coulomb drag, as represented in 
Fig. [l] is the phenomenon where a voltage V2 applied to a 
two-dimensional (2D) conducting layer (called the active 
layer) generates both a current I2 on that plane and a volt- 
age V\ in another layer, parallel to the first, and located 
at a distance d 1 . This effect occurs because electrons in 
the active layer, under the presence of an electric field £?2, 
"drag" the electrons in the other layer through their mu- 
tual Coulomb interaction. The current densities in each 
layer, ji and j'2, are related to electric fields in each layer 
through the conductivity tensor, according to the general 
relation: 



(i) 



where an is the conductivity of each isolated layer (d — >• 
00) and o~rj is the so-called trans-conductivity. Notice that 
in a drag experiment no current flows in the non-active 
sheet, that is, ji = 0, and hence we can express the elec- 
tric fields Ei and E2 in terms of the current 22 in the 
active plane alone. This allow us to define the quantities 
of experimental interest, namely the drag resistivity 
given by 

WV 1 E, a D 
I2L j 2 deter 
and the longitudinal resistivity pL, reading 

E2 an 

32 



where deta = an 022 — cr% is the determinant of conduc- 
tivity tensor a (W is the width and L is the length of the 
device) . 
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Fig. 1: Generic representation of a Coulomb drag experiment. 

Using the above equations we can eliminate the trans- 
conductivity in terms of pjj and rewrite the relationship 
between the drag and the longitudinal resistivities as: 



PL 



Ipl/Pi 



where pi = 1/au. Eq. Q has two extreme limits: 



(4) 
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PL ~ ^ \Jp2lp\pD , PD > VP1P2/2 . ^ 
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Hence, in the limit where the drag resistivity is small 
compared with the resistivity of the two layers (pr> <C 
y/piP2/2) the longitudinal resistivity is dominated by the 
resistivity of the active layer. This is the usual case of the 
2D electron gases (2DEG), found in semiconducting het- 
erostructures 2-4 , where the drag resistance is a small 
effect. 

In what follows, we will argue that for the case of an 
ultra-clean graphene double layer the opposite regime can 
be reached, namely, one can have pD ^> \Jpip2/2 with 
Pi <^ P2 so that the longitudinal resistivity, p^, is propor- 
tional to the drag resistivity, but enhanced by a factor pro- 
portional to \J pijp\ > 1, so that p L > p D > y /p 1 p 2 /2. 
This regime can be reached in ultra-clean graphene by 
tuning the bottom and top gates in the device (see Fig. 
[I]) such that the Fermi energy of the non-active layer, 
epi = vpkpi (where vp ~ 10 6 m/s is the Fermi veloc- 
ity), is large, but the Fermi energy of the active layer, 
6p2 = ^f^F2, is small and close to the neutrality point 
(that is, one requires that kpi ^> kp2)- This is only possi- 
ble in graphene because the resistivity in each individual 
graphene layer is a monotonically decreasing function of 
the density, |n e |, in each plane (kpi — \J^\n e \) 5 . In 
ultra-clean graphene, the resistivity is a very sharp func- 
tion of the density and therefore the longitudinal resistiv- 
ity can be enhanced by orders of magnitude by the drag 
effect. In fact, we show that the drag resistance diverges 



at low densities as po 



-0.8 



(n e — » 0) and therefore the 



longitudinal resistance can become arbitrarily large. We 
stress that the behavior for pp, as computed using our 
model, is better described by the relation p£> ~ n e a<y7le \ 
where the exponent a(n e ) is a function of n e and shows a 
crossover from a ~ 2 to a ~ 1 as n e decreases. 

The mechanism described here can produce a huge 
enhancement of the on-and-off ratio for current flow in 
graphene, without the need for the opening of a gap in 
the spectrum, solving the famous bottleneck for the use of 
graphene devices in high end electronic applications. 

While we wait for the drag data in ultra-clean samples 
[6] , we check our model against the existing data for drag 
resistivity measurements in graphene on Si02 u\. As it is 
well known, Si02 is a dirty substrate, and at low densities 
electron- hole puddles are formed [M9] , greatly limiting the 
mobility in these devices. We show that our model applies 
to this conventional case describing the data extremely 
well. 

We notice that the theoretical literature on Coulomb 
drag in graphene is scarce [TolITT], and it was shown [7] 



that the theoretical approaches used so far 10|[TT 



are un- 
able to describe the experimental data. Under a number of 
simplifying assumptions, it was shown lTj that electrons 
in graphene, when described by the massless Dirac equa- 
tion, should have zero drag-resistivity. According to that 
analysis 11 , trigonal warping corrections would be nec- 



dependence of the scattering time is taken into account 
and dynamic effects in Coulomb screening are correctly 
included. 

Theory of Coulomb drag and experimental re- 
sults. — The drag resistivity can be obtained from the 
solution of Boltzmann's kinetic equation 



12 



16 . In this 



approach it is assumed that the main scattering mecha- 
nism within a graphene layer is electron-impurity scatter- 
ing 17 18 and that the electronic density is outside the 



range where one finds electron- hole puddles 8,9|. In dirty 
substrates such as SiC>2 this can happen at densities of 
10 12 cm -2 , however, in cleaner substrates such as Boron- 
Nitride [l9] this only happens at extremely low densities 
of the order of 10 10 cm -2 20 . Without loss of general- 
ity we also assume that the graphene layers are electron 
doped and, as explained above, the Fermi energy in the 
two graphene layers is such that e F 2 < cfi- We also make 
use of the full dynamical screening between the layers, 
which takes into account intra- as well as inter-layer inter- 
actions. However, in calculating the resistivity we only 
take into account intra-band interactions between elec- 
trons belonging to each of the sheets (the validity of this 
assumption is achieved by keeping the electronic density 
large enough [II] ) . Finally, one key point of our approach 
is taking into account the full momentum dependence of 
the electronic scattering time tj, , originated from electron- 
impurity scattering. It is well known that is roughly 
proportional to the square root of the electronic density, 
that is, we have tj, = rofc, where To is a constant com- 
puted elsewhere [17], and which drops out at the end of 
the calculation. This assumption is in agreement with the 
experimental data [5] and is essential for the accurate de- 
scription of Coulomb drag in graphene. 

Within these assumptions, Boltzmann's kinetic equa- 
tion suffices for the description of the drag resistivity. The 
final result for the latter quantity (and the central result 
of this paper; see derivation ahead) is 



1 y/eFltF2 2xY7 7 m i \ 

Pd = a c F{k F 2,k F i,T,d c ,a g ) 

go z°7rkB-L 

= -poJ 7 (k F2 , k F1 ,T, d c ,a g ) . 



essary to explain a finite drag resistivity in graphene. We 
show here that such is not the case when the momentum 



(6) 



The reader is referred to the derivation of Eq. ( |24| ) for 
the definition of the several quantities in Eq. (|6|. Using 
Eq. © we are able to describe quantitatively the drag 
resistivity measured in graphene at different temperatures 
[7]. Our results are summarized in Fig. [5] In that fig- 
ure the symbols are the experimental data and the dashed 
lines are the calculated values of pr>- The drag resistivity 
is represented as function of the gate voltage. The calcu- 
lated curves have been horizontally shifted by Vd = — 2 V, 
since the graphene in the device had its neutrality point 
at V g = Vd, whereas the calculated curves assume neutral 
graphene when V g = 0. The agreement is quantitative up 
to V g = 4; below we discuss the origin of the deviations for 
V g < 4 V, which are associated with the electron-puddle 
formation in the device. To understand the experimental 
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a=1 .7, a c =0.21 , a b =1 .2, d =14 nm, e =6 




Fig. 2: Calculation of the drag resistivity as function of the 
gate-voltage and comparison with the experimental measure- 
ments 7 , using a c — 0.21, db — 280 nm, 6b — 3.8, dt — 14 
nm, and e t = 6. Note that except for a c , the parameters db, 
dt, 6b, and et are a characteristic of the apparatus used in the 
experiment, and therefore fixed numbers. The symbols are the 
experimental data and the dashed lines are the calculated val- 
ues of po- The inset represents the electronic density, n\ and 
712, in the two sheets, as obtained from the solution of Eqs. |7| 
and (J8|. We note that m and 712 are not independent of each 
other. 



conditions associated with the measurements reproduced 
in Fig. [2] the reader is referred to Ref. [7]. 

The device has two different dielectrics: Si02 (relative 
permittivity 65 = 3.8), with thickness d b = 280 nm, and 
AI2O3 (relative permittivity e t = 6 [21]), with thickness 
d t = 14 nm (we note that other authors report a value of 
e Ai 2 o 3 =9-1 [22] ). The latter oxide is sandwiched between 
the two graphene layers. As detailed in Ref. 7 , the den- 
sities in the two layers are not independent of each other. 
The electronic densities n 2 and n\ follow from [7]: 



eV g = d b (n 2 +ni) + v F h^Jim 2 , (7) 



d t m 




(8) 



The numerical solution of the last two equations gives the 
electronic density in the two graphene sheets for any value 
of V g (see inset in Fig. |2j. 

In our theory, the only free parameter is the intrinsic 
dielectric of constant graphene or, equivalently, the inter- 
action parameter: 



2.2. 



(9) 



We note that there is some controversy on the value of 
otg [23]. When graphene is immersed in a dielectric of 
relative permittivity e r , the interaction parameter be- 



are three different interaction parameters: two associ- 
ated with electron-electron interactions within the two 
graphene planes (a t and a&) and one associated with the 
inter-plane interaction (a c ). The constants a t and a b play 
their role in the calculation of the dielectric function of the 
coupled layers. Explicitly, we have: 

ott = 2a g /(l + e A i 2 o 3 ) , (10) 

Qlc = O^/eAlaOsJ (11) 

a b = 2a g /(e Ah o 3 + e Si o 2 ) • (12) 

Taking a c as reference, we get a t ~ 1.7a c and a b ~ 1.2a c ; 
these two numbers are fixed in our model. Therefore, the 
only free parameter is a g , or in alternative a c . We have 
chosen a c as our fitting parameter and the best fit was 
obtained for a c = 0.21. This latter value implies: 



a„ ~ 1.3 . 



(13) 



a number close to the one given by Eq. 

The vertical dashed line in Fig. [2] sets the limit of va- 
lidity of our theoretical model in what concerns its ap- 
plication to graphene on dirty substrates. Indeed, at that 
point the experimental curves change their curvature from 
negative to positive, indicating a change of regime. Below 
V g = 4 V, the theoretical results start to deviate from 
the measured data. This happens because the electronic 
density in the bottom layer approaches the regime where 
electron-hole puddles control the electronic transport in 
graphene 7 . In fact, looking at the green curve of Fig. 
3. a in Ref. 7 , we see that during the transition from the 
regime of electron-electron drag to that of electron-hole 
drag the behavior of po is reminiscent of that reported 
for Hall-measurements in graphene close to the neutral- 
ity point [5]. In both cases, the elementary theory, which 
ignores the effect of electron- hole puddles, predicts diverg- 
ing drag resistivity and Hall coefficient. We note that our 
computed curves diverge as the electronic density is re- 
duced; the divergence is enhanced at high temperatures. 
In the absence of electron-hole puddles pD can become 
large, giving rise to the switching effect we have described 
in the introduction. 

Derivation of the drag resistivity formula. — The 

calculation of the current density 22 (the driven current) 
and the electric field E\ (the drag field) requires the solu- 
tion of Boltzmann's equation, under the assumption that 
the current density in sheet 1 is zero, j\ =0. The driven 
current 22 reads 



32 = E 2 go7rv F ron 2 . 



(14) 



Following a standard procedure 
is given by: 



14 



16 , the drag field E\ 



£1 



1 e 1 g s g v f dki 
g tttii J 4tt 2 



hi cos 0fc x 



df ku 
dt 



(15) 



comes a r 



j/e r . In the experimental setup 7 there 



where dfk liC /dt] ee is the collision integral due to inter- 
plane Coulomb interactions, e\ is the charge of the carriers 
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in plane 1, go = 2e 2 //i, rt\ is the electronic density in 
layer 1, and g s and g v are the spin and valley degeneracy, 
respectively. We stress that unlike previous works [To] 
we take fully into account the momentum dependence of 
the relaxation time. It was noted before that ignoring 
this momentum dependence leads to incorrect results [l6] 
(for the 2DEG, however, the relaxation time can be taken 
momentum independent). 

The calculation of the collision integral requires the use 
of the Coulomb interaction between electrons in different 
layers. The latter is given by 



H 



12 



1 

4A 



£ £ U dc (q)fa 
a k 1 a a ki+qa'bl l2 pb k2 - q p> , 



'(ki,q)fpp>(k2,q) x 
(16) 



where the creation operators and b k2(X refer to elec- 
trons in layers 1 and 2, respectively, A is the area of the 
system, and 



eie 2 ^ 
2e c e q 



-qd c 



V c (q)e 



-qd c 



(17) 



where d c is the inter-layer distance. The transferred mo- 
mentum between the layers is q — k\ — k 2l a, a', /?,/?' = 
v, c refer to the valence (y or -1) and conduction (c or +1) 
band electrons, the momentum of electrons in plane 1 is 
denoted by k\ and for electrons in plane 2 by k 2 \ the mo- 
mentum sums h\ and k 2 contain spin summations as well. 
The chiral nature of the electron wave function is encoded 
in the form factors: 



foLOL'{k\,q) = 1 
fpp'(k 2 ,q) = 1 

where = arctan^//^). 



3 ) 



(18) 
(19) 



Following a standard approach 14 16 , 24 , we have 
dfci,. dfk u 



J 

I 4^1 



4tt 2 
hduj 



hi cos 6 kl 



dt 

\U(q : uj)\ 2 



smh z (huj/2k B T) 



_ _ E 7rev F r 
2 2hk B T X 

e - 

P lc (q,uj)P 2c (q,uj) , 



The difference between our calculation and that of Tse et 

(i) 

at. 10 is precisely in the form of the G k q . The dynamic 
dielectric function e(qr, uj) for the two-layer system is given 
by [I6j[25]: 

e(q, uj) = e 2 (q, u)ei(q, uj) - Uj c (q)Vi(q, uj)V 2 (q, uj) , (23) 

where €j{q 1 uj) = 1 — Vj(q)Vj(q,uj) and Vj(q,uj) are the 
dielectric and polarization functions of each of the indi- 
vidual layers [26]. Notice that a& and a t enter in e^q^uo) 
and e t (q,uj), respectively. Once the functions Pj C (q,uj) 
are computed, a relatively simple expression follows for 
the drag resistivity: 



Pd 



1 \f^F\^F2 

go 2 5 7rk B T 



a 2 c J r (k F2 ,k FU T 1 d c1 a g ) , (24) 



where U(q,uj) = Ud c (q)/e(q,uj) is the dynamic Coulomb 
interaction, evaluated using the random phase approxima- 
tion for the dynamic dielectric function e(q,uj). We note 
that ro is not the relaxation time (see above), and it can- 
cels out when the ratio E\/j 2 is computed. We also have 



dki 



G k j ijf kj , q 5(huj + huj kjjq ), (21) 



4^2 kj,q 



%-q/2 



where 6f kjtq = f kj+q / 2 - f kj - q / 2 , 

%+ q /2 (j = 1,2), and the function G^ q , with j = 1,2, 
is defined as 



hu k 



\kj - q/2\cosO kj _ q/2 



2 + 2*i , j 

\k 3 +q/2\cos6 k3+q/2 \] \f cc (k J ,q)\ 2 .(22) 



which leads to Eq. (|6| and a c is the interaction parameter 
of graphene considering the dielectric in between the two 
graphene layers made out of AI2O3. The dimensionless 
function J r {k F2 , kpi,T, d c , a g ) is defined as: 



F(k F2 ,k F1 ,T, d c ,a g ) 



JL r Xm r y 

m=l J ° Jo 



where 
F(x,y) 



dyF(x, y) , 

(25) 



x 9 (x 2 -y 2 ) $ 1 (x,y)$ 2 {x,y) 
\<x,y)\ 2 sinh 2 ^h^k j^y 



-2d c x\/kFik£ 



(26) 

and x and y are dimensionless variables defined as q — 
\Jk F \k F2 x and uj = v F \/k F ik F2 y, respectively. The three 
integration regions are: < x\ < 6, < yi < min(x, — x + 
b); b/2 < x 2 < (a + 6)/2, max(-x + 6, x - b) < y 2 < 
min(x, — x+a); a/2 < £3 < 00, max(— x+a, x—b) < ys < q 
where a = 2^/k F \/k F2 and b = 2^/k F2 /k F \. The function 
e(x, y) is defined as: 



e(x,y) 



(x^/x 2 - y 2 + P 12 )(x\/x 2 -y 2j r P 2 \) 



Pi2P2ie 



(20) where 



4a. s 



-2d c \fkF\k,Fix 



.a 



(27) 



■ly^^y), (28) 



and a s equals or a t , depending on the graphene layer. 
Below we outline the calculation of the functions $ s (x, y). 

In Fig. [3] we plot the quantity poF(x,y) for particular 
values of k F i (as determined by the gate voltage), the ef- 
fective interaction parameter a c , the distance between the 
two layers, and the temperature. In that figure the three 
panels, from top to bottom, represent the contribution 
of the three integration domains defined by x m and y m . 
At low temperatures and high back-gate voltage, only the 
contribution of region m = 1 is significant; on the other 
hand, at high temperature and low densities all the three 
regions (m = 1, 2, 3) give a significant contribution to p£>, 
as shown with the example in Fig. [3J 



p-4 



Coulomb Drag and High Resistivity Behavior in Double Layer Graphene 




£> F(x,y) 




Fig. 3: Plot of the integration kernel poF(x,y) in Eq. (25) for 



the parameters V g = 3 V, d — 14 nm, a g — 0.21, and T — 300 
K. The three figures correspond to three different regions of 
the (g, ft) plane. From top to bottom the regions are ordered 
as m =1 1 m =2, and m =3 (see text). 



Particular limits of pp. — Let us now consider the 
case where a&, a c , and a t are all equal, as is the electronic 
density in the two sheets. Taking the limit of very low 
densities, we see that pr> can reach values of the order 
of 10 3 ft, as in the case represented in Fig. [Z] It is the 
possibility of having large values of po that suggests the 
mechanism for turning on-and- o/f the electric current in a 
double-layer graphene-based transistor. We also note that 
in the density range n e ~ 3 x 10 12 cm -2 to n e ~ 0.5 x 10 12 
cm -2 , pD is better described by pD oc rie a ^ ne \ where 
the exponent a(n e ) shows a crossover from a(n e ) ~ 2 to 
a(n e ) ~ 1 as the electronic density is reduced. 

Finally, when e F ^> fceT and d c k F ^> 1, the asymptotic 
regime for p£> reads 



Pd oc- 



1 {k B Tf 



go v F 



(29) 



The asymptotic formula (29) has a different dependence on 



d and n e from that found in Ref. [To] , due to the inclusion 
of the momentum dependence in the scattering time; in 
Ref. [lO] the integration over q has a kernel of the form 
q 3 / sinh (qd c ), whereas our kernel, in the same limit, has 
the form q b / sinh 2 (g<i c ). Thus we obtain pr> oc d 6 and Tse 
et al. have obtained pD oc d A . The difference by a factor 
of two in the power of the momentum q comes from: (i) 
r q oc q; (ii) r q appears twice in the right-hand side of Eq. 



'q_f 



15k. It is worth stressing that Eq. rt29fc does not apply 



to the conditions of the experiment we have described in 
Fig. [2j and the experimental conditions where it would 
apply can hardly be reached. Also, the values of pn of 
an experiment done in the regime of validity of Eq. ( 29 ) 
would be difficult to measure. 

Technical details: Calculation of the func- 
tions <frj(x,y) at zero temperature. — The functions 
$>j(x,y) - j = 1,2 - depend on the integration limits x\ 
and yi - i = 1,2,3; we also define £1 = uo/vp. We write 
$j(q,n,k Fj ) = $f(q,Sl,k Fj ) + $j(q,n,k Fj ), where 

/oo 
dcj) sinh 2 (j)fj [(q cosh 2 ± ft) /2] . 
-co 

(30) 

At zero temperature we have: 

^[(gcosh^ift)/^] = 0[2k Fj - (q cosh 2 <j>± ft)] . (31) 

In terms of the functions ®j~(q, £l,kpj) and $j(q, ft, k F j) : 
the function $>j(q,Q,k F j) reads: 

1. region 1: < q < k F j and < ft < g, in which 
case we have <£>j = — $~ (we have omitted the 
arguments of the functions &j(q, ft, kpj) for the sake 
of simplicity); 
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2. region 2: k F j < q < 2kpj and < Q < 2kpj — <?, in 



which case we have 



3 3 



3. region 3: k F j < q < 2kpj and 2kpj — q < ^ < in 



which case we have 



-*7 



4. region 4: q > 2kpj and q — 2kpj < £1 < <?, in which 



case we have <E>j 



-*7- 



In the regions where &f is finite we have 



<s>+( q ,n,k Fj ) 

<f>-(q,fl,k Fj ) 



-^min 

+ - sinh(20 min ) , (32) 



7 max 

+ - sinh(2</> max ) , (33) 



where (/> m ax = arccosh [(2k pj + fy/ol an d ^min = 
arccosh [(2k pj — &)/q]- The four regions above have their 



ft = 2k F - q 



O = 2k F -q 




Fig. 5: Domains where the functions <&~t (q,Q,kFj) (left) and 
<E>~ (q,Q,kFj) (right) are finite. In the triangular area, both 
functions are finite. In the rectangular area only one of the 
functions is finite. 

graphical representation in Fig. [5j In the triangular area 
of the diagram [5| the polarization reads 



Vj(q,n) 



2k 



Fj 



= <&,(<?, O). (34) 



In the rectangular area of the same diagram, the polariza- 



tion is obtained by adding to Eq. (34) the term 

— *— —=l — e,-(g,n), 

where Oj(q, 0) = O a (q, 0) + ©b(g, 0) and 

9 a (<j,0) = (2^-0)^9-2^+0 



(35) 



1/9 + 2fc F j - O , 



6 (g,0) 



-2q arctan A 





n 


- 2k Fj 




n 


+ 2k Fj 



(36) 
(37) 



We have used these results for the polarization in the cal- 
culation of the dynamic dielectric function. 

Conclusions. — We have given a quantitative the- 
ory of Coulomb drag in graphene. We have shown that 
for ultra-clean graphene, unlike the case of a conventional 
2DEG, can be tuned to a region of diverging drag and 



longitudinal resistivity allowing the use of double layer 
structures for device applications where the on-and-off ra- 
tio for the current flow has to be large. We also show 
that our theory explains quantitatively the experimental 
data for the drag resistivity in dirty devices away from 
the electron-hole puddle region. As this region shrinks in 
cleaner devices, the validity of our theory extends. 

Supplementary Information. — In Fig. [6] we plot 
the dependence of the drag resistivity as function of the 
electronic density, n e , the interlayer distance, d ci and the 
temperature, T. The density range scanned in this figure 
is different from that given in Fig. [4j In the three panels 
of Fig. [6] the black dashed lines are fits to a power law 
of the form po ~ x ax , where x represents one of three 
parameters n e , d ci and T. From these fits we have found 
that pD follows, in the regime of parameters of Fig. [6| 
roughly the behavior 



Pd 



(Tk B ) 2 

v F h 2 nldl 



(38) 



It is also clear from Fig. [6] that Eq. (38) holds only ap- 
proximately. Also, comparing the dependence of po on 
n e , as given in the regime of parameters of Fig. [4j with 
that of Eq. (38), obtained from Fig. |6j we see that the ex- 
ponent a n& changes with density, being smaller at smaller 
density values. 



k F =0.08 -- 0.41 nm 




1 2 3 4 5 6' 6 8 10 12 14 20 40 60 80 100 

n (x10 12 /cm 2 ) d c (nm) T (K) 



Fig. 6: Dependence of the drag resistivity as function of the 
electronic density, n e , the interlayer distance, d c , and the tem- 
perature, T. The used value of a is a = 0.21 
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